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Abstract. Fast algorithms for aritlimetic on real or complex polynomials are well- 
known and have proven to be not only asymptotically efficient but also very practi- 
cal. Based on Fast Fourier Transform, they for instance multiply two polynomials 
of degree up to n or multi-evaluate one at n points simultaneously within quasi- 
linear time 0(n ■ polylogn). An extension to (and in fact the mere definition 
of) polynomials over fields R and C to the skew-fieXA H of quaternions is promis- 
ing but still missing. The present work proposes three approaches which in the 
commutative case coincide but for EI turn out to differ, each one satisfying some 
desirable properties while lacking others. For each notion, we devise algorithms 
for according arithmetic; these are quasi-optimal in that their running times match 
lower complexity bounds up to polylogarithmic factors. 



1 Motivation 

Nearly 40 years after CooLEY and TuKEY |4|, their Fast Fourier Transform (FFT) 
has provided numerous applications, among them 

- fast multiplication of polynomials 

Given the coefficients of p,q C[X], n :— deg(p) + deg{q); 
determine the coefficients of p ■ q. 
which, based on FFT, can be performed in (n • log n) and 

- their multi-evaluation 

Given the coefficients ofp e C[X], deg(p) < n, and xi, . . . , a;„ G C; 
determine the values p{xi), . . . ,p{x„). 
allowing algorithmic solution within 0{n ■ log n). 

Observe in both cases the significant improvement over naive (n^ ) approaches. These 
two examples illustrate a larger class of operations called Fast Polynomial Arithmetic 
I1I14I with, again, a vast number of applications |7J- For instance, Gerasoulis em- 
ployed fast polynomial arithmetic to drastically accelerate A^-Body Simulations in 2D 
[8 1, and Pan, Reif, and Tate did so in 3D 1 11 1. Since systems with up to = 10'"' 
objects arise quite frequently when simulating biochemical processes, the theoretical 
benefit of asymptotic growth 0{N ■ poly log A^) over 0{N'^) pays off in practice as well. 

Technically speaking in order to calculate, for each of the N particles, the total force 
it experiences due to the — 1 others, Gerasoulis identifies the plane with C; 
he thus turns Coulomb's potential into a rational complex function which, by means of 
fast polynomial multiplication and multi-evaluation, can be handled efficiently. 1.1 1..13J 
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on the other hand exploit fast multi-evaluation of polynomials to approximate the total 
forces in M.^. Whether 3D forces can be obtained exactly within subquadratic time is still 
an open question. One promising approach proceeds by identifying, similarly to fsl, 
with (a subspace of) Hamilton's four-dimensional algebra of Quaternions H and there 
applying fast polynomial arithmetic of some kind or another In fact the mere notion 
of a polynomial becomes ambiguous when passing from fields IK = M and IK = C to 
the skew-field K = H. We consider three common approaches to define polynomials 
(Section|2j and, for each induced kind of quaternion polynomials, present quasi-optimal 
algorithms supporting according arithmetic operations (Section|3}- 

2 Quaternions 

The algebra H of quaternions was discovered in 1843 by W.R. HAMILTON in an attempt 
to extend multiplication of 'vectors' from = C to R'^. In fact, H is a/oMr-dimensional 
real vector space whose canonical basis l,i,j,k satisfies the non-commutative multi- 
plicative rule 

— j"^ ^ k'^ — ijk = —1, ij = —ji — k + cyclic interchange (1) 

which, by means of associative and distribute laws, is extended to arbitrary quaternions. 
IH is easily verified to form a skew-field, that is, any non-zero element a possesses 
a unique two-sided multiplicative inverse a^^. In fact it holds a^^ = a/|ap where 
a :— Re(a) — i Imi(a) — j Imj(a) — k Imfe(a) is the analogue of complex conjugation 
and \a\ :— \/a ■ a — \/a ■ a E M+ the norm satisfying \a ■ b\ = \a\ ■ \b\. The center 
of HI is K; in other words: real numbers and only they multiplicatively commute with 
any quaternion. For further details, please refer to the excellent' CHAPTER 7 of |5|. 
Theorem 17.32 in p| determines the (multiplicative algebraic) complexity of quater- 
nion multiplication; |2 1 does so similarly for quaternion inversion and division. However 
rather than on single quaternions, our focus shall lie on asymptotics w.r.t. n, the quater- 
nion polynomials' degree, tending to infinity. 

It is well-known that commutativity has to be abandoned in order to turn R"* into 
some sort of a field; in fact, Frobenius' Theorem states that HI is the only associative 
division algebra beyond ^ C. On the other hand to the author's best knowledge, all 
notions of polynomials either require the ground ring 31 to satisfy commutativity or — 
such as skew polynomial rings, see P.262, CHAPTER 16 of [10 1 — they lack evaluation 
homomorphisms. The latter means that any polynomial p = p{X) e ^[X] should 
naturally induce a mapping p : ^ Jl, x i-^ p{x) such that for all a, a; G 3?: 

X(x) — a{x) — a, fr~q{x) ~ p{x) ■ q{x) , and p + q{x) = p{x) + q{x) . 

The distant goal is to find a notion of quaternion polynomials which naturally gen- 
eralizes from real or complex ones and supports efficient arithmetic by means of, say, 
quasi-linear time algorithms. Our contribution considers three such definitions for K[X] 
which, in case K is an infinite field, are equivalent to the usual notion. In case K = H 

' wrongly condemned in CHAPTER XXI, P.245 of QH. . . 
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however they disagree and give rise to different arithmetic operations. We focus on Mul- 
tiplication and Multi-Evaluation and present in Section|3] for each of the three notions, 
according quasi-optimal algorithms. 

2.1 Polynomials as Ring of Mappings 

The idea pursued in this subsection is that the following objects should be considered 
polynomials: 

- the identity mapping X := id : IK ^ K, a; i— > x, 

- any constant mapping a : IK ^ K, x ^ a for a e IK 

- the sum of two polynomials and 

- the product of two polynomials. 

Formally, let the set IK"^ of mappings / : K ^ K inherit the ring structure of K by 
defining pointwise f + g : x ^ f {x)+ g{x), f -g : x ^ f{x)-g{x).T\\tntmbtd 
K into this ring by identifying a G IK with the constant mapping IK 9 a; a £ IK. 

Definition 1. IKi [X] is the smallest subring o/K^ containing X and the constant map- 
pings IK. For instance, 

ai+X ■a2- X ■ X ■a'i + a4^- X ■ X ■ X -az eKi[X] , ai, ... £K fixed. (2) 

Ki [X] is closed not only under addition and multiplication but also under composition, 
i.e., f + g,f-g,fog & ^i[X] for f,gE Ki[X]. Since, in the commutative case, any 
such polynomial can be brought to the form 

n e N, Of e IK , (3) 

Definition [3 there obviously coincides with the classical notion of polynomial rings 
M.[X] and C[X]. For the skew-field K = H of quaternions, the structure of IHIi[X] is 
not so clear at first sight; 

- a-X ^ X- a unless a e R i.e., the form (|3} in general cannot be attained any more. 

- Uniqueness becomes an issue, since 

X -X -i-X-i + i-X-X-i-X - i-X-i-X-X - X-i-X-X-i (4) 
vanishes identically |5 TOP OF P.201]; 

in particular, a polynomial can have many more roots than its 'degree' suggests. 

- The fundamental theorem of algebra is violated as well: i- X ~ X -i + l has no root 
inH0P.2O5]. 

- Lagrange-style polynomials Pm to pairwise distinct points xq, . . . , Xn-i G H, e.g., 

( n — 1 s ^ n — 1 \ —1 TL—1 

Y[{X~Xi)y(l[{x^-xe)] or [{x^- x,)-^ ■ {X - x^)'^ 

both interpolate Pm{xm) = 1, Pm{xi) — 0, m ^ £ but obviously lack uniqueness. 

- There is no polynomial division with remainder; e.g. X-i-X mod = ??? 
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On the other hand we present in Subsection l3.2l algorithms for addition, muhiphca- 
tion, and multi-evaluation of this kind of quaternion polynomials of degree n in time 
0{n^ ■ polylogn). Since it turns out that generic p £ iii[X] have roughly free co- 
efficients, the running time is thus quasi-optimal. Finally, a fast randomized zero-tester 
for expressions like (|2j and comes out easily. 

2.2 Polynomials as Sequence of Coefficients 

Since the above Definition ^ thus does not allow for quaternion polynomial arithmetic 
as fast as quasi-linear time, the present subsection proposes another approach. The idea 
is to identify polynomials with their coefficients. Recall that for p = X]"=o^ ae^^ and 
q = J2"=o^ biX^ over a commutative field K, the finite sequence of coefficients c = 
(q) e k* of p • g is given in terms of a = (a^) £ IK* and b ~ (bi) <E K* by the 
convolution product 



with the implicit agreement that a/ = for £ > n and be ^ for £ > m. 

Definition 2. K2[X] is the set K* of finite sequences of quaternions, equipped with 
componentwise addition and convolution product according to 0. Let X denote the 
special sequence (0, 1, 0, . . . , 0) € K*. 

It is easy to see that this turns K2 [X] into a ring which, in case of fields K of character- 
istic zero, again coincides with the usual ring of polynomials K[X]. Here the classical 
results assert that arithmetic operations + and * can be performed within time 0(n) and 
0(n • logn), respectively. In Subsection 13. II we show that the same is possible in the 
non-commutative ring ]Hl2[X]. Dealing with n coefficients, this is trivially quasi-optimal. 

Unfortunately fast arithmetic for H2 [X] does not include multi-evaluation, simply 
because evaluation (substituting X for some a; G H) makes no sense here: One might be 
tempted to identify a e H* with the formal expression a^X^ and b with b^X^, 
but then c a * b does not agree with 



{E<^iX')-{EbeX') = EiELo<^fX'-b,^fX'-' ^ Y.iC,X' 



because of non-commutativity. 

The next subsection considers expressions of the form ^ a^X^ as further notion 
of quaternion polynomials. These lack closure under multiplication; on the other hand, 
there, multi-evaluation does make sense and turns out to have classical complexity {n ■ 



2.3 One-sided Polynomials 

Roughly speaking, one aims at a subclass of ]HIi[X] where polynomials have only 0(n) 
rather than 0{n^) coefficients and thus give a chance for operations with quasi-linear 
complexity. 



c = a * 6, ce ^Y.t=o°-t ■ b£-t, £ ^ 0, ...,n+m-l 



(5) 
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Definition 3. Let X : IK ^ K denote the identity mapping and consider this class of 
mappings on K; M.z[X] := { YTt^o a^X^ : n G No, &£ G K} C K^. 
The degree of p e ^^[X] is deg(p) = max£, deg(O) := -1. 

Again this coincides for fields K of characteristic zero with the usual notions. For the 
skew-field of quaternions, the restriction compared to Q applies that all coefficients ai 
must be on the left of powers X^. Unfortunately, this prevents H3 [X] from being closed 
under multiplication; fortunately, IEIl3[X] has the following other nice properties: 

- being a real vector space; - allows fast multi-evaluation; 

- supports interpolation; - a fundamental theorem of algebra holds; 

- polynomials satisfy uniqueness. Formally: 

Lemma 4. Consider p :— X)"=o^ ciiX^, ai G H. 

a) Suppose p{x) ~ Ofor all x G H. Then a£ = Ofor all £. 

b) Nevertheless even p 7^ may have an infinite ( and in particular unbounded in terms 
ofp's degree) number of roots. 

c) If ai ^ Ofor some £ > I, then p has at least one root. 

Proof, a) Follows from Lemma Eb) by choosing n > deg(p) and pairwise distinct 
xq, . . . , Xn-1 G K since then, no three are automorphically equivalent. 

b) All quaternions x — i(3 + jj + kS with (3,'y,S G M and [3'^ + j"^ + 5'^ = 1 are easily 
verified zeros of p :— X^ + 1. 

c) Cf. P. 205 in |3 or see, e.g., 0. 

Interpolation is the question of existence and uniqueness, given xq, . . . ,Xn-i and 
yo, . . . , Un-i G K, of a polynomial p G IK[X] with degree at most n — 1 satisfying 
p{xi) = ye for all ^ = 0, . . . , n — 1. In the commutative case, both is asserted for pair- 
wise distinct X(. Over quaternions, this condition does not suffice neither for uniqueness 
(LemmallJ)) nor for existence: 

Example 5. No p = aX^ + bX + c e M^lX] satisfies p{i) = = p{j), p{k) = 1. 

It turns out that here an additional condition has to be imposed which, in the commu- 
tative case, holds trivially for distinct xi, namely being automorphically inequivalent. 

Definition 6. Call a, 6 G H automorphically equivalent iff a = u-h ■ for some 
non-zero u G H, that is, iff Re(a) = Re(6) A |Im(a)| = |Im(6)| where 
Im(a) :— ilmi(a) + j lmj(a) + fclmfc(a). 

This obviously is an equivalence relation (reflexivity, symmetry, transitivity). The name 
comes from the fact that mappings x ^ u ■ x ■ are exactly the M-algebra automor- 
phisms of H; cf. L5, BOTTOM OF P.215]. The central result of |9J now says: 

Lemma 7. For ccq, . . . , ^n-i G H, the following are equivalent 

a) To any yg, . . . , j/n-i G U, there exists p G IHl3[X] q/deg(p) < n such that p{xe) — 
ye,£^0,...,n-l. 

b) Whenever p = J2l=o ^^^^ '^"'^ 1 — X]"=o ^^^^ satisfy p(xe) — q{xi) for £ = 
0, . . . , n — 1, follows ae = bg. 
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c) The Quaternion Vandermonde Matrix V := (a;™)j!,,„=o, « invertible. 

d) Its Double Determinant \\V\\ does not vanish. 

e) The X£ are pairwise distinct and no three of them are automorphically equivalent. 

Concluding this subsection, ]Hl3[X] has (unfortunately apart from closure under mul- 
tiplication) several nice structural properties. In 13.31 we will furthermore show that it 
supports multi-evaluation in time 0(n • log^ n). More generally, our algorithm applies 
to polynomials H|[X] { X;"Jo^ • • &f : n e No, 5^ G H} 

with coefficients to both sides of each monomial X^. This generalized notion has the 
advantage of yielding not only an K-vector space but a two-sided H-vector space. 

3 Algorithms 

3.1 Convolution of Quaternion Sequences 

Beginning with the simplest case of H2 [X]: 

Let n e N. Given a = (ao, ai, . . . , a„_i) G H"- and b = (&o, 61, • ■ • , ^m-i) G H™, 
one can compute their convolution according to (|5} from 16 real convolutions^ and 12 
additions of real sequences within time 0(n • logn). Indeed write componentwise 

a = Re(a)+i Imi(a)+j' Im^ (a)+fc Imfc(a), 6 = Re(6)+i Imi(6)+j Imj(6)+fc Imfc(fa) 

and exploit M-bilinearity of quaternion convolution. 

3.2 Ring of Quaternion Mappings 

The central point of this subsection is the identification of Hi [X] with the four-fold 
Cartesian product of four-variate real polynomials Y\ R[Xo, Xi, X2, X3]. Formally 
consider, for / : H — > H, the quadruple / of four-variate real functions defined by 

/o(Xo,..,X3):=Ro(p(Xo + iXi+jX2 + fcX3)) /i(Xo,..,X3):=Imi(p(Xo + iXi+iX2+fcX3)) 
/2(Xo,..,X3):=Imj(p(Xo + iXi+iX2+fcX3)) /a (X,,. .., X3):=Im;.(p(Xo + iXi +iX2 +fcX3)) 

and multiplication among such mappings /, 5 : M"* ^ M"' given pointwise by 

(/o,/l,/2,/3) • {90,91,92,93) ■= (6^) 

(/oSo-/iSi-/2(;2-/3S3, /03l+/lS0+/2S3-/3S2, f 092+ 12 90+ fa 91 - f 193 . /0S3+/3S0+/lS2-/23l) 

4 

In that way, calculations in Hi [X] can obviously be as well performed in Yl ^[Xq, .., X^]. 
This allows for application of classical algorithms for multivariate polynomials over 
commutative fields. But before, we need a notion of degree on Hi[X]: 

Definition 8. For a commutative multi-variate polynomial, let deg denotes its total de- 
gree; e.g., deg(x^i/^) = 5, deg(O) = ~oo. The degree deg(q) of a quaternion polyno- 
mial q € Hi[X] is half the total degree of the real four-variate polynomial /q + . . . + 
/I with fo, ■ ■ ■ , fs according to (|6}. 

^ In fact, 4 complex convolutions suffice; but asymptotically, that gains nothing. 
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Rather than the total degree, one might as well have considered the maximum one 
deg{x^y^) :— 3 since, for 4 variables, they differ by at most a constant factor How- 
ever we shall later exploit the equality deg{p ■ q) = deg(p) + deg{q) valid for the first 
whereas for the latter in general only the inequality deg{p-q) < deg{p) + deg{q) holds. 
In fact, this nice property carries over to the degree of quaternion polynomials: 

Lemma 9. The degree deg(p) of p G Hi [X] is always integral. Furthermore it holds 
deg{p ■ q) = deg(p) + deg{q). 

Now recall the following classical results on four-variate polynomials: 

Lemma 10. a) Given (the coefficients of) p,q £ C[Xq, . . . , X3], the (coefficients of 
the) product p ■ q can be computed in time 0{n* ■ logn) where n := deg{p ■ q) = 
deg(p) + deg{q). 

b) Given p £ C[Xo, . . . , X^\ of degree n, one can compute within 0(n^ • logn) steps 
the coefficients of p(T ■ {Xq, X^)^ + 2/) G C[Xo, X3], that is, perform on p 
an affine variable substitution given by T G C*^^ and y G C^. 

c) A given polynomial p G C[Xo, . . . , X^] of degree n :— deg(p) can be evaluated on 
all n"^ points of a 4-dimensional complex grid G := Aq x Ai x A2 x A3 such that 
Ag (Z C \Ae\ — n, within time 0(rt^ • log^ n). 

d) The same holds for the regular affine image G' = T ■ G + y of such a grid, i.e., 

G' ^ {T-a; + y : a; = (xo,. .. ,a;3)t G G}, T e C^""* regular, y e . 

e) Let p G C[Xo, . . . ,^3] be non-zero, n > deg(p). Fix arbitrary A C_ (C of size 
\A\ > 2n. Then, for {xq, . . . , x^) G A* chosen uniformly at random, the probability 
ofp{xo, . . . , ^3) — is strictly less than ^. 

Proof, a) Reduction to the univariate case by means of Kronecker's embedding: 
cf. Equation (8.3) on r62 of | 1 | for m := 4; dealing with the complex field C 
rather than an arbitrary ring Jl of coefficients, the loglog-factor may be omitted. 

b) Folklore. A proof had to be removed from the final version due to space limitations. 

c) Cf. Equation (8.5) and the one below on P.63 of 1 1 1 for m := 4, c := n. 

d) follows from b). It is not known whether multi-evaluation is feasible on arbitrarily 
placed points within time 0(?i^ • polylog7i). 

e) Cf. Subsection 12.1 in|14|. 

One could of course identify in a similar way complex univariate polynomials p G C[Z] 
with tuples Po,Pi G M.[X, Y] of real bivariate polynomials. However the thus obtained 
running times of 0(n^ • poly logn) thus obtained for C[Z] are strikingly suboptimal, 
basically because not every tuple of real bivariate polynomials corresponds to a com- 
plex univariate polynomial. For instance, z i— > Re(z) is well-known not only to be no 
complex polynomial but to even violate Riemann-Jacoby's equations of complex dif- 
ferentiability. Surprisingly for quaternion polynomials, the situation is very different: 

Lemma 11. Re{X) = i(X - iXi - jXj - kXk) G Hi[X]. More generally, every 
quadruple of real four-variate polynomials corresponds to a quaternion polynomial. 
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The generic quaternion polynomial of degree n thus has 0{n'^) free coefficients. Lem- 
mas ^|and[^together yield 

Theorem 12. a) Multiplication of two quaternion polynomials p, q G ]HIi[X] is possi- 
ble in time C)(n^ • logn) where n :— deg(p • q) — deg(p) + deg(q). 

b) Multi-evaluation of p at a;o, ■ • ■ ,Xn-i S H can be done within 0(r7,'* • log^ ri), 
n := deg(p). 

c) Within the same time, multi-evaluation is even feasible at as many as points x, 
provided they lie on a (possibly affinely transformed) n^-grid G. 

The above complexities are optimal up to the (poly-)logarithmic factor. 

Theorem[21presumes the polynomial(s) to be given as (coefficients of four) real four- 
variate polynomials. But how fast can one convert input in more practical format like (|2j 
or to that form? By means of fast multiplication of several polynomials, this can be 
done efficiently as well: 

Theorem 13. a) The (ordered!) product Pt of ra quaternion polynomials pi G 
Hi [X], each given as quadruple of real four-variate polynomials, can be computed 
within 0{ri^ ■ logn ■ logm) where n — Aeg{pi) denotes the result's degree. 

b) An algebraic expression E over quaternions, i.e., composed from +,—,-, constants 
a € H, and the quaternion variable X — but without powers like nor brackets! 
— can be converted into the quadruple of real four-variate polynomials according 
to ^ within time 0{N^ ■ log^ N) where N = \E\ denotes the input string's length. 

The above conversion yields a deterministic 0{N^ ■ log^ A^)-test for deciding whether a 
given quaternion expression like (0} represents the zero polynomial. When satisfied with 
a randomized test, the same can be achieved much faster: 
TheoremlOlfcontinued) 

c) Given £ > and an expression E of length N = \E\, composed from "+", 

" ■ ", constants a € H, the quaternion variable X, and possibly brackets " ( ", " ) "; 
then one can test with one-sided error probability at most e whether E represents 
the zero-polynomial within time 0{N ■ log j-). 

Proof a) Standard divide-and-conquer w.r.t. m similar to Corollary 2. 15 in Js). 

b) Lacking brackets, the input string E necessarily has the form 

E = El it E2 lb . . . zb E]\4 

where Eg describes a product Pi of quaternion constants (degree 0) and the inde- 
terminate X (degree 1). Since obviously deg(P£) < Ni := \Ei\, its real four- 
variate representation is obtainable within G{Nf ■ log^ Ni) steps. Doing so for all 
£ = 1, . . . , 7\f leads to running time 0{N^ ■ log^iV) as < ^■ 

W.l.o.g. let deg(Pi) < deg(F2) < ■•• < deg (Pa/)- Adding up the just ob- 
tained four-variate representations in this increasing order takes additional time 
0{Nf + + ... + Nfj) < 0(iV4). 

c) By virtue of standard amplification it suffices to deal with the case e = i. The 
algorithm considers any set yl C M of size |^| > 2N. It chooses xq, xi, X2, X3 G A 
uniformly and independently at random; and then evaluates the input expression E 
by substituting X := xq + ixi + jx2 + kx^. If the result is zero, the algorithm 
reports zero, otherwise non-zero. 
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The running time for evaluation is obviously linear in iV. Moreover, only one- 

sided errors occur So suppose E represents non-zero p e Hi[X]. Then obviously 
deg(p) < N and at least one of the four real four-variate polynomials pQ, . . . 
according to (|6j is non-zero as well. By virtue of LemmallOb). this will be witnessed 
by {xq, . . . , X3) — i.e., p{xo + ixi + jx2 + kx^) 7^ — with probability at least 



3.3 Multi-Evaluating Two-Sided Polynomials 

Consider an expression of the form p{X) = aeX^be, ai, bg S H. Expanding 

ae = Ke{ai) + ilmi(a^) + jlmj{ai) + klmk{ae) and similarly for b^, one obtains, 
by virtue of distributive laws and since whole R commutes with X^, that it suffices to 
multi-evaluate expressions of the form 

q{x) = Ero«^^'' a,eR{l) (7) 

since p{X) can be obtained from 16 of them, each multiplied both from left and right 
with some basis element 1, i, j, k. Now with real a£, multi-evaluation of is of course 
trivial on xq, . . . , Xn-i G C; but we want X£ to be arbitrary quaternions! Fortunately, 
the latter can efficiently be reduced to the first. 

To this end, consider mappings (pu : M ~^ M, x t-^ u ■ x ■ u^^ with u e EI of norm 
|m| = 1. It is well-known |5 PP.214-216] that, identifying H with M"^, ipu describes a 
rotation, i.e., tpu G SO(K"'). Furthermore, restricted to the set 

ImH {a; e H : Re(x) = 0} ^ 

of purely imaginary quaternions, exhausts whole SO(R'^) as u runs through all unit 
quaternions; this is called HAMILTON'S Theorem. Finally, ipu is an (and in fact, again, 
the most general) M-algebra automorphism, i.e., satisfies for a G M and x,y E M: 

ipuia) = a, (pu{x + y) = ipuix) + (fiuiy), (pu{x ■ y) = (pu{x) ■ (puiy) ■ 



Lemma 14. Forv, w G ImH, |?;| = 1 = \w\, letu := {v+w) /\v+'w\; then (pu(v) = w. 
In particular for X G H\R, v := Im(x)/| Im(a;)|, w := i, it holds q{x) = ■q{y)-u 
where y := u ■ x ■ u^^ G K + = C. 

Our algorithm evaluates q G M[X] simultaneously at xi, .., a;„ G H as follows; 

- For all a;£ G R + iK, let ui 1; 

- for each xi ^ M, compute (in constant time) ug according to Lemma [141 

- Perform in linear time the transformation yi := ug ■ xg ■ uj^ . 

- Use classical techniques to multi-evaluate q at yi, ...,?;„ G C within 0(n • log^ n). 

- Re-transform the values q{yg) to q{xe) — uj^ ■ q{yg) ■ ug. 



This proves the claimed running time of 0(n • log^ n). 



□ 
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4 Conclusion 

We proposed three generalizations for the notion ' polynomiaV from fields R and C to 
the skew-field H of quaternions and analyzed their respective properties. For each no- 
tion, we then investigated (where applicable) on the algebraic complexity of operations 
multiplication and multi-evaluation on polynomials in terms of their degree. The upper 
bounds attained by our respective algorithms match (usually trivial) lower bounds up to 
polylogarithmic factors. 

However since each of the above notions lacks one (e.g., closure under multiplica- 
tion) or another (e.g., quasi-linear complexity) desirable property, a satisfactory defini- 
tion for quaternion polynomials is still missing. Here comes another one, generalizing 
the representation of complex polynomials in terms of their roots: 

Ki[X] := {ao • (X - ai) • (X- as)- • -(X -a„) : n G No,a^ e H} (8) 

So what is the complexity for multi-evaluation in IHI4 [X] ? 

In view of the planar A^-body problem, Gerasoulis' major break-through was fast 
multi-evaluation of complex rational functions 



N 



(9) 



1=1 



for given oi, . . . , aM G C at given xi, . . . , xjv G C; cf. also COROLLARY 7 in fTS^. 
Our techniques from Subsection 13. 31 vield the same for a:;^ G H and ai G K. Thus 
the crucial question remains whether (|9j also allows multi-evaluation in sub-quadratic 
time for both and xi being quaternions. But what is a rational quaternion function, 
anyway? We do not even know what a quaternion polynomial is! Observe that, lacking 
commutativity. 



X-a X-b X~a X~b ' ' ' ' X-a X-b 

cannot be collected into one single fraction, in spite of the common denominator. 

Acknowledgments: The idea to use quaternions for iV-Body Simulation was suggested 
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A Postponed Proofs 

Here, we collect some proofs which, in the printed version, had to be quelled due to page 
constraints. 

Lemma 15. Let /, g G R[^i, . . . , Xd\ denote d-variate real polynomials. Then, 

deg{f+g^) = max{deg(/2),deg(g2)} = 2 max{deg /, deg^} . 
In particular, the total degree of + is even. 

Proof. The second equation is trivial, regarding that the total degree satisfies deg(/^) = 
2dcg(/); similarly for the inequality deg(/^ + g'^) < max{dcg(/^), dcg(g^)}. One 
thus has to show that, although in + g^ certain terms of coinciding maximum total 
degree might indeed cancel, the above inequality is in fact an equality. This is where 
the real ground field comes into play: choose in / and g respective terms M and N 
of coinciding maximum total degree, that is, M — a ■ X"^^ ■ ■ ■ X^'' and N = b ■ 
X^' ■ ■ ■ X^-" with ^ = deg / = deg g = ^ nt and a,b ^ Q. Then both = 
a2.j^2mi . . .^2mrf ^ 52.^2«i . . . x^"" have total degree cqual to deg(/2) = 

deg{g^). Furthermore, their respective occurrences in + g^ cannot cancel because 
and are stricdy positive; hence deg(/^ + g^) > deg(M^) = deg(-/V^). □ 

Proof of Lemma |5| Integrality of the degree is covered by applying Lemma [Tsl above 
inductively to [f^ f?), {fg + fl /|), and (/^ + ff + /|, /|). 

For the second claim, straight-forward calculation confirms the Four Squares Theorem^ 
for real numbers, here applied to the case of real polynomials /o, . . . , /a, go, . . . j 53^ 

(/o + fl + fl + fl) ■ {gl + ~9l + ~9l + 33) = 

= (/offo - - 7232 - /353)^ + (/o5i + /i5o + /253 

+ (/off2 + /2go + /3gi - fihf + Uoh + /350 + /152 
Now observe that the total degree of the left hand side of (I10> is 

deg {{fl + fl + fl + fl) ■ ((.go' + ~9l + ~9l + .93)) 

= deg(/2 + fl + fl + fl) + deg(g2 + ~gl + gl + ~gl) 

which, by Definition|8] agrees with 2dcg(p) + 2deg(g) forp, g e IHIi[X] according 
to (|6}. At the same time, in view of Equation (6^), the total degree of (fTol's right hand 
side is nothing but 2 deg(p • g). □ 

Proof of Lemma nn 

• Straight forward calculation verifies Re(X) — j{X — iXi — jXj — kXk) which 
obviously belongs to Hi [X] . Thus, the quadruple of four-variate real polynomials 
{Xq, 0, 0, 0) e n K[Xo, . . . , X3] does correspond to a quaternion polynomial. 

^ discovered by Euler in 1748 



- /352)' (10) 

- /251)' 
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• Similarly, (Xi, 0,0,0) G Jl^I^I^o, • • • , ^3] corresponds to Iin,(X) = Re{-iX) e 
Hi [X]; same for {X2, 0, 0, 0) and (X3, 0, 0, 0). 

• For any real constant a, {a, 0, 0, 0) G Yi^ ^[Xq, . . . , X3] corresponds to a £ H C 

Mi[X]. 

• Let (/, 0, 0, 0) and (5, 0, 0, 0) correspond to quaternion polynomials p, g e Hi [X], 
respectively. Then (/ + g, 0, 0, 0) corresponds to p + q E Mi [X]; and, 

(/ -3,0,0,0) ^'i^ (/, 0,0,0) -(5, 0,0,0) 
corresponds to p ■ q eMi [X] . 

Since IR.[Xo; . . . , X3] is the smallest set containing real constants, the generators Xq, Xi, 
. . ., X3, and being closed under addition and multiplication, the above considerations 
imply that, for any / E R[Xq, . . . , (/, 0, 0, 0) corresponds to some p e Hi [X]. 

• Suppose (/, 0,0,0) corresponds to p G Hi[X]. Then (0,/,0, 0) corresponds to 
-ip G mi[X]; analogously for (0, 0, /, 0) and (0, 0, 0, /). 

• Let (/o, 0,0,0) correspond to po e Mi[X], (0, /i, 0, 0) topi, (0, 0, /2, 0) top2, and 
(0, 0, 0, /s) to p3. Then (/o, /i, /s, h) corresponds to pq + . . . + ps E Mi[X]. □ 



Proof of Lemma [T3I Observe that u ^ = tt/|ttp = — u since u E ImH and |u| = 1. 
Thus 

{v + w)v{v + w) + v^w + wv^ + wvw 

= [;;T^^ " 2 + 2{v,w) 

because |t;| = 1 = \w\ by presumption. As G ImH, furthermore = — 1 = -up 
and thus 

—v — 2w + (2(— w, w)w + v) 
2 + 2{v, w) 

The R-algebra homomorphism property ensures that q((^„(x)) = ipu{q{x)) for any 
polynomial q with real coefficients and x G H. In particular, M-linearity yields 

y ^ ipuix) = i^„(Re(a;) + I Im(a;)| • u) = Re(a;) + | Im(a;)|i G M + iR 

□ 



